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1. Introduction 



This work is a continuation of the study of semileptonic B decays by the HPQCD collaboration 
||l|], which was announced last year. Since then, we have included all 1/M corrections for the 
temporal and spatial parts of the heavy-light electroweak currents |Vb| and |Vyt|. Previously, only 
the corrections that are lowest order in 1/M were included. Also, last year our data was on one 
ensemble, now we have more ensembles, listed in and hence more full QCD data points in 
addition to partially quenched results. Moreover, we have started implementing and examining the 
effects of staggered chiral perturbation theory fits using formulas by Aubin and Bernard [^. 



2. Current Matching 

We constructed the heavy-light currents in terms of lattice operators so that they accurately 
represent the continuum heavy-light currents up to the desired order. We picked a process where the 
operator to be matched is relevant. The corresponding scattering amplitude was matched between 
lattice and continuum through one loop in perturbation theory. The result is independent of the 
process picked to do the matching. Our process is a heavy quark scattering into a light quark 
through the electroweak current. The corresponding Feynman diagram is calculated through one 
loop in continuum perturbation theory; then the result is expanded in powers of 1/M. The terms 
in the expansion are identified as the continuum analogs of the necessary lattice operators. These 
lattice currents are: 



Temporal Currents: 



jS'\x) =q{x)nQ{x), 



Spatial Currents: 



where = 7^ or 757^, and T, = 1 or 75. 

Because of the symmetries of our action, axial vector and vector currents give the same results. 
We used this as a further check in our matching calculations. 

The continuum currents are a linear combination of the lattice currents through 0{as/M): 
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{V,),ont = (l + «5P0)(4*") + (l + «5Pl)(4'^'™') + «5P2(4'^'™'), 

mean, = (1 + «5P0)(/r) + (1 + «5Pl ) (4'^'™') + (^SPl^J^^'"") + «5P3(4'^) + «5P4(4''), 

where 7q ''"'^ = 7q ' — as^mJ^^^ . The subtracted currents /j,'' "'^ are more physical, as they have 
power law contributions subtracted out through 0{as/aM). 

Last year, our calculation was correct through C?(o;.v): it only included the first term, (1 + 
(XsPo) {j'^^)- We now take into account all the 1 /M terms in the above expression. To this end, we 
calculated both the nonperturbative matrix elements (7) and the perturbative coefficients p,. The 
perturbative calculations are described in 

The difference between matching at 0{as/M) and at 0{as) is found to be very small. This 
difference is plotted for one of the currents, 7q^^ in Fig.l. 
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Figure 1: The size of 1 /M current correction is shown. The more physical subtracted currents have an 
even smaller effect than the naive ones, with the central value of the lai^gest /jj^^ around 1 percent. 



3. Form Factors 



The semileptonic form factors are defined by the following relation with the hadronic matrix 
element of electroweak currents between a B meson and a pion: 
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{K{pn)\V^\B{ps)) = U{q^) 



where v'' = ^, p^^ = p'^ - {pj,-v)v^, q>'=p1^-p'^. 
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Figure 2: Linear extrapolations of f± to the physical chiral limit at several fixed values of E,i. 
The three-point correlator 

d'\p,jB,t,TB) = Y^Y^{^n{0)J^{lt)^l{y,T^^^^^ 

z y 

where O's are interpolating operators for the mesons, is fit using Bayesian fitting algorithms to the 
following form to obtain Aqq: 

d'\p.,PB,t,TB)^ I I V^^'(^-'). 

Aqo is directly related to form factors /j| and in the following way: 

Aoo(V°) 
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^Tt bS Pit 

where Aqo includes the effects of matching, and are the groundstate amplitudes from the B 
meson and the pion 2 — pnt correlators respectively. 



4. Procedure for Approaching the Chiral Limit 
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Figure 3: The form factors f+[q^) and fo{q^), calculated using the same ensemble as the previous year, 
including the l/M current corrections are shown and compared with the results from last year. The curve 
shown is a BK [^] fit to last year's data. Results from other ensembles and more sophisticated chiral extrap- 
olations have not been included yet in this plot. 



One needs to extrapolate the form factors obtained as described in the previous section to the 
physical chiral limit. We first interpolate /|| and f± to the same values of Ej^. After interpolating, 
we do chiral extrapolation at several fixed Ej^ values to get the form factors at the physical point. 
A plot of f± showing the interpolated points at fixed Ej^, and linear extrapolation to the physical 
limit can be seen in fig. 2. So far, only simple linear fits were used in our results. We have started 
looking into doing more complicated chiral fits using 
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Using the chirally extrapolated /ji and /j^, we obtain f+{q^) and fo{q^)- The form factors 
f+{q^) and fo{q^), calculated using the same ensemble as the previous year, including the l/M 
current corrections are shown in Fig. 3 and compared with the results from last year. The curve 
shown is a BK [^] fit to last year's data. Our results are in the region of q^ larger than l5GeV. 
We perform fits of our results to model ansatze for the q^ dependence in order to obtain the form 
factors at smaller q^. A new approach "Moving NRQCD" shows promise to directly calculate 
the form factors even at small q^. 

5. Summary 

Last year, unquenched studies of semileptonic B decays using NRQCD heavy and improved 
staggered (AsqTad) light quark actions were reported. The main change since then is that we have 
included all l/M corrections to currents through 0(a^,aas, a^/aM, a^Agc/j/M), instead of only 
0{as). We also have more ensembles, instead of just one ensemble with dynamical quark mass 
am = 0.01. We have started coding and examining the effects of staggered chiral perturbation 
theory. The results shown in this article, however, used only linear fits. We also have new data on 
finer lattice spacing that needs to be fit and included. Once our analysis is complete, and we have 
the corresponding f+{q^) and fo{q^), this will enable us to extract the CKM matrix element \Vi,h\ 
using the new experimental B— > nlv branching fraction results, for instance due to Belle, BaBar, 
CLEO. 
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